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Preface 



There is a deep analogy between the physics of crystaUine sohds and the behavior 
of superfluids, dating back to pioneering work of Philhp Anderson, Paul Martin and 
others. The stiffness to shear deformations in a periodic crystal resembles the super- 
fluid density that controls the behavior of supercurrcnts in neutral superfluids such as 
He^. Dislocations in solids have a close analogy with quantized vortices in superflu- 
ids. Remarkable recent experiments on the way rod-shaped bacteria elongate their cell 
walls have focused attention on the dynamics and interactions of point-like dislocation 
defects in partially ordered cylindrical crystalline monolayers. In these lectures, we 
review the physics of superfluid helium films on cylinders and discuss how confine- 
ment in one direction affects vortex interactions with supercurrcnts. Although there 
are similarities with the way dislocations respond to strains on cylinders, important 
differences emerge, due to the vector nature of the topological charges characterizing 
the dislocations. 
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Introduction 



Defects of various types have long been known to play a key role in condensed matter 
physics and materials science. Point-like vacancy and interstitial defects in crystalline 
solids facilitate atomic diffusion and give rise to color centers (|T|). Vortex lines in type 
II superconductors, singularities that concentrate magnetic field lines, move in the 
presence of an external electrical current ([HE]), causing dissipation unless the vortices 
are strongly pinned, an especially difficult problem for practical applications of high 
temperature superconductors, due to thermal fluctuations Soft metals bend irre- 
versibly when a dilute concentration of translational line defects called dislocations (see 
below) move in response to an external stress ([S]). The Bronze Age began when early 
metallurgists added small concentrations of tin to copper to create a stronger and more 
effective material - the tin impurities pin dislocations in place, thus preventing plastic 
deformation. The Iron Age was the result of making iron even stronger than bronze 
by supplementing the addition of impurities (typically carbon) with work hardening 
at the blacksmith's forge, which entangles the dislocation lines (jS]). Finally, melting of 
two-dimensional crystalline solids with increasing temperature can be driven by the 
serial unbinding of thermally activated point-like dislocations and then point-like rota- 
tional defects called disclinations, with an intervening hexatic phase (j7j)- The detailed 
theory found early confirmation with experiments on well-equilibrated 2d assemblies of 
colloidal particles (|5]), and most recently in a remarkable decade- long sequence of ex- 
periments on super- paramagnetic particles interacting with a repulsive 1 /r^ potential 
by the group of Georg Maret at the University of Konstanz (|9|). 

Remarkably, dislocation dynamics also plays a significant role in the two-dimensional 
covalently bonded macromolecules that comprise the cell walls of ^ 1 micron-diameter 
rod-shaped bacteria. A fragment of such a cell wall is shown in Fig. II. lb . consisting of 
a regular meshwork of glycan {i.e., sugar) strands of alternating moieties of NAM (N- 
acetylmuranimic acid) and NAG (N-acetylglucosamine) , with these chains connected 
by short peptide {i.e., amino acid) crosslinks (fTOj) . This peptidoglycan meshwork ranges 
from a single layer in some thin-walled gram-negative bacteria to of order 8-15 layers 
in the thicker cell walls of gram positive organisms. Because of its nonzero shear mod- 
ulus, the cell wall maintains the cylindrical shape of the bacteria, as well as resisting 
anywhere from 1-20 atmospheres of outward osmotic pressure. If the cell wall is gently 
eliminated (and rupture is avoided), rod-shaped bacteria become spherical, due to the 
fluid character of the confinement by the remaining lipid bilayers. While carrying out 
their important functions, cell walls must maintain their integrity as bacteria such as 
E. coll elongate to twice their initial size approximately every 20-30 minutes, followed 
by cell division. It is as if a pressurized dirigible aircraft had to be systematically 
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Fig. 1.1 (a) Peptidoglycan layer of a bacterial cell wall. Glycan strands of alternating NAM 
and NAG sugars are connected by short peptide linkages, (b) An edge dislocation arising 
from the insertion of a extra half-row of unit cells like that in (a) ; as the dislocation moves 
upwards, it mediates easy extension of a glycan free end and is characterized by the Burgers 
vector 6, the amount by which a circuit that would close in a perfect crystal fails to close due 
to the defect. 

remodeled to gradually become twice its initial size while still in flight! 

It turns out that point-like edge dislocation defects provide a natural way for two- 
dimensional partially ordered bacterial cell walls to remodel themselves. As illustrated 
schematically in Fig. II. lb . an edge dislocation can be regarded as an extra half-row 
of unit cells ([S] [TT|) . Dislocations are characterized by a topological charge called the 
Burgers vector b, given by the amount the contour integral of the displacement field 
u{x, y)around the defects fails to close, 



du ~ b. 



(1.1) 



If the extra row of unit cells terminating in dislocations like those in Fig. 11.21 could 
be systematically extended upwards by adding new material, the crystal would widen 
along the axis of the cylinder. The cylinder will elongate at constant radius as a sub- 
set of "activated" dislocations circumnavigates via this climbing motion and extends 
glyan strands. This defect mechanism for bacterial elongation, with the activated de- 
fects called "glycan extension centers", was proposed by Burman and Park in 1984 
(fT2l [T3| , approximately 50 years after the importance of dislocations was recognized in 
materials science ([T4) [TSl [T6|) . From a materials science perspective, it is at first hard to 
see how dislocation dynamics could ever be relevant in networks of the strong covalent 
bonds that characterize a bacterial cell wall - dislocations were invented to understand 
plastic deformations of metals, where atoms are connected by much weaker metallic 
bonds. How can the very strong bonds in Fig. II. lb break and reform to remodel the 
cell wall at physiological temperatures? Biology accomplishes this seemingly impossi- 
ble task because it has evolved a specialized suite of enzymatic molecular machines 
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Fig. 1.2 Schematic of dislocations on the cylindrical portion of a bacterial cell wall. Circum- 
ferentially rotating activated dislocations propelled by attached strand elongation machinery 
are indicated by arrows. Inactive dislocations are indicated by asterisks. In reality, the in- 
activated dislocations are likely to greatly outnumber the activated ones. Here, Ro is the 
constant radius the cylindrical cell wall of the bacteria, and L{t) is its length, which grows 
as dislocation motion causes the bacterium to elongate. 

with names such as MreB, RodA, PBPlb, PBP2, MreC and MreD There is a 

rough analogy between this collection of enzymes and the enzymes that allow DNA 
replication every time a cell divides (|12p. Instead of facilitating the complementary 
base pairing necessary to complete a single strand of DNA, these enzymes remodel 
cell walls by breaking peptide crosslinks, extending glycan strands and then adding 
new amino acid cross-bridges. Interestingly, only ^ 20-30 out of the approximately 
10,000 glycan strand ends in a typical bacterial cell wall ([T^ need to be activated 
to allow a cell to double its length in about 30 minutes ((T^)- Pairs of broken strand 
ends are associated with the end of a row of unit cells like those in Fig. II. lb . and thus 
qualify as a dislocation, according to the criterion illustrated in Fig. II. lb . The suite 
of remarkable enzymes that creates an activated dislocation can only lengthen of the 
cell wall - new enzymes are needed to form the septum that partitions a growing cell 
into two daughters and completes a cell division ([T7| . 

In these lecture notes, we will not describe in detail the application of dislocation 
theory to this fascinating biophysics problem. A detailed exposition can be found 
in Refs. (|18l [T9l [20)) . Instead, we focus on the simpler, but mathematically closely- 
related problem of point-like superfluid vortices disrupting order in the condensate 
wave function of superfluid helium films on cylinders. Superfluid helium has a rich and 
illustrious history, enriched by ideas from such giants of theoretical physics as Fritz 
London, Lars Onsager and Richard Feynman (pil . The deep analogy between 

the physics of crystalline solids and the behavior of superfluids and superconductors, 
first fruitfully exploited by scientists such as P. W. Anderson and P. C. Martin ([24l 
I25p . is summarized for two dimensional systems in Table I. Here, u{'r) is the same 
phonon displacement field measuring deviations from a perfect crystalline lattice that 
appears in Eq. (|l.ip . and 6{r) describes the phase of the superfluid order parameter. 
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Its relation of its gradient to the superfiuid velocity in terms of Planck's constant H 
and the mass m of a helium atom will be discussed further below. The energy of both 
systems is quadratic in the derivatives of these fields, with coefficients or "stiffnesses" 
given by elastic constants and the superfiuid density respectively. The last row of the 
table is the quantization condition on point-like defects, given by a Burgers vector b 
equal to a lattice displacement for dislocations in crystalline solids, and by an integer 
s for vortices in superfiuid helium films. 



Table I: Analogy between the elastic deformations of two-dimensional crystalline solids 
and low energy excitations of superfiuid He^ films. 





Crystalline Solid 


Superfiuid Helium 


Dynamical field 


strain matrix 

"u(^) = 1 [dxUj{r) + djUi{r)] 


superfiuid velocity 


Energy 


E^hJd'r [2^iul + Xul,] 


u = \psjd\ \^ve{f)[ 


Stiffness 


^ and A (Lame' coefficients) 


Ps ^ {Ps ~ superfiuid density) 


Defect contour integral 


f du{r) = b 


§de{r) = 2tts, s = 0,±1,±2,... 



In what follows, we describe the phenomena associated with the right-hand column 
of Table I, eventually showing how the excitation of quantized vortex pairs triggers the 
decay of super currents on cylinders, thus extending the earlier work of Ambegaokar et 
al. on planar helium films (|26p . In Sec. II, we review the physics of inviscid superfiuid 
helium films. As we shall see, low energy excitations of these very thin films, which can 
be less than a 100 atomic layers thick, resemble gravity waves in a classical fiuid of finite 
depth, except that the dominant restoring force is the van der Waals interaction with 
the substrate. These van der Waals waves, also called "third sound", were predicted 
and observed experimentally by K. R. Atkins and coworkers approximately 50 years 
ago ([27l [28]) . Third sound fails to propagate above temperatures exceeding a few 
degrees above absolute zero ([29l [30| . In the limit of long- wavelength third sound, 
corresponding to approximately uniform supercurrents, this failure can be attributed 
to the excitation of quantized vortex pairs with equal and opposite vorticity. In Sees. 
Ill, and IV we show how the energetics of superfiuid vortices can be mapped onto 
two-dimensional electrostatics embedded in a cylinder, with a crossover from a 2d 
logarithmic potential at short distances to a linear confinement when distances exceed 
the cylinder circumference. After discussing the nuclcation of vortex pairs on a cylinder 
in the presence of a mixture of longitudinal and azimuthal supercurrents, in Sec. V 
we briefiy compare and contrast the energetics of the vector dislocation charges for 
two-dimensional crystals wrapped around a cylinder. Charge-neutral dislocation pairs 
with Burger's vectors at a generic angle to the cylinder axis behave much like their 
counterparts in scalar 2d electrostatics. However, for the special case of Burgers vectors 
aligned with the cylinder axis, relevant to elongation mechanism for bacteria shown 
in Fig. II. 2[ the interactions fall off exponentially at large distance (|19l f^O)) . due to a 
remarkable screening effect associated with vector charges, related to the physics of 
grain boundaries (jS]). 
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Third Sound: van der Waals waves 
in thin helium films 



Below the lambda temperature Tx = 2.17°Kelvin, liquids composed of He*, the bo- 
son isotope of helium, flow without resistance between the narrowest parallel plates, 
through the smallest capillary tubes and through tightly packed powders like jewelers 
rouge. Unless extreme care is taken, these fluids can easily defy gravity and climb 
out of their confining vessels (PT|) . Remarkably, this superfluid phase (also denoted 
He II), fails to crystallize and remains a liquid down to absolute zero at atmospheric 
pressure. At temperatures below Tx, the dynamical properties of the superfluid phase 
can be understood in terms of a "two-fluid model" , consisting of a normal fluid density 
Pn{T), with a conventional shear viscosity, and a superfluid density ps{T), for which 
this viscosity is entirely absent (|5T|) . He II exploits its remarkable ability to flow with- 
out resistance by wetting nearly all substrates, with the exception of materials with 
very weak van der Waals interactions such as Cesium ([5^ Thus, as illustrated 
schematically in Fig. 12. 1[ a pool of bulk superfluid liquid helium in a closed chamber 
(in equilibrium with its non-superfluid vapor phase) will coat all walls with a superfluid 
film, including the cylindrical pedestal in the center. 

If we focus for the moment on the thin film on the circular top of the pedestal at 
height H above the bulk fluid, the normal fraction will be completely pinned by its 
viscosity, due to the close proximity of the substrate. The normal fluid density will, 
in any case, typically be much less than the superfluid density in the low temperature 
limit considered here. The superfluid velocity Vs{f,t) within this film will then obey 
the inviscid Navier-Stokes (or Euler) equation, 

-^ + ivs- V)i7, = Vp + f,,t, (2.1) 

at ps 

where p{r,t) is the fluid pressure and fext{r) = — V$(r) is a constant external body 
force due to a potential energy function $(r) acting on every parcel of fluid. Incom- 
pressibility implies that V • iT, = 0. What is the thickness d of this film? In the limit 
H ^ d indicated in Fig. 12.11 the potential energy $o exactly at the top of the film at 
the point B, relative to the point A, arises from a combination of the gravitational 
potential and the energy that binds helium atoms to the substrate. 



<^,o{H,d)^mgH-^, 



(2.2) 
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Fig. 2.1 (a) A pool of bulk superfluid helium (He II) in closed chamber with cylindrical 
symmetry coats all walls with a superfluid helium film, including the cylindrical pillar in the 
center. The pillar projects a height H above the surface of the bulk liquid, and the film has 
thickness d. (b) Van der Waals interaction energy (Ifvdwiz) for a helium atom at height z 
above a substrate. This diverging energy leads to very high pressures for small z, causing a 
dead layer of inert non-superfluid helium atoms to form at distances less than do from the 
substrate. 

where m is the mass of a helium atom. The second term arises from the van 
der Waals attraction of a helium atom at the film-vapor interface to the substrate, 
integrated over the 3-dimensional half-space occupied by the pillar. The constant /? 
depends on the polarizability of the material that forms the pillar pip . Note that the 
zero of the gravitational energy is set at the height of the bulk liquid and the zero of 
the interaction with the substrate is for a film thickness d — oo. In equilibrium, the 
value of $0 at a point B on the surface of the bulk fiuid well away from the pillar 
must be exactly the same as at point A, which leads to a prediction for the thickness 
of the film, 

d{H) ^ w 4 X 10^^(1 cmjUfl'^. (2.3) 

where the numerical estimate follows from the value of /3 for a glass substrate P7|) . The 
prefactors for other substrates are similar. For a macroscopic pillar height H 1 cm, 
the film will be only d ~ 200 helium layers thick! 

As indicated in Fig. 12.1b . 4>vdw {z) , the van der Waals contribution to the potential 
energy of a liquid parcel at height z above the substrate has a short distance cutoff 
at a distance do- For d < do, the substrate forces are so large that the helium atoms 
solidify into a "dead layer", in which superfluidity is absent. A superfluid film thus 
only exists for a range of heights z with do < z < d. Upon inserting gravitational and 
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Fig. 2.2 A van der Waals wave at the liquid- vapor interface of a superfluid helium film with 
equilibrium thickness d. This excitation resembles a gravity wave, except that the dominant 
restoring force is the van der Waals interaction with the substrate. 

substrate interaction energies into Eq. (|2.1[) , we find that the velocity field in the film 
on top of the pillar at height z obeys 



dt 



1 



(2.4) 



Consider the relative magnitudes of the two downward body forces, fg = —gz and 
fvdw = —{3l3/z'^)z on the right hand side of Eq. (|2.4p . Upon taking z ^ d{H) = 
, we havefydw/fg ~ '^P/d'^g = iH/d. Thus, although the gravitational 
energy helps set the height of the film, it is utterly negligible compared to the van der 
Waals forces when H ^ d. 

Consider the dynamics of the ripple-likc third sound excitation at the liquid- vapor 
interface, as indicated schematically in Fig. 12.21 This wave involves both a distortion 
of the interface, and a nonzero velocity field below the surface. In Fig. 12. 2[ both the 
interface height h{x,z,t) and the velocity field are independent of y. As mentioned 
above, third sound is similar to a gravity wave in a classical fluid (|34]), except that 
the restoring force of gravity replaced by interactions with the substrate. Like grav- 
ity waves, the motion is potential flow, i.e., Vs{x,z,t) = "S/ 4>{x, z,t). However, unlike 
classical fluids, for important quantum mechanical reasons, the potential function is 
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related to the phase 9{x, z, t) of the wave function that describes the macroscopic Bose- 
Einstein condensation that underUes the dissipationless superfluid ([U |3T]) . Indeed, at 
sufficiently low temperatures, we have 

Vsix,z,t) ^ —e{x,z,t) = V(t){x,z,t), O<0<27r, (2.5) 
m 

Where m is the mass of a helium atom, and the restriction on the phase angle is 
natural because its relation to the condensate wave function is approximately "00 ('^) ~ 
ypje'*^''-' (|55]). With the neglect of gravitational forces, Eq. (|2.4p can thus be written 

where we have neglected the nonlinear advective term (iT, • 'V)vg. Let us assume the 
excitation in Fig. 12.21 has period r and wavelength A. A criterion for the neglect of the 
nonlinear coupling follows if we first parametrize the interface height as 

h{x, t)^d + C{x, t) = d + Coe'fc"-*"* (2.7) 

and then note that, dvs/dt ^ Co/t^ whereas {vs ■ V)vs ^ Co/''' ^- Our neglect of 
advection is thus justified for Co <SC A, i.e. for excitations with wavelengths much 



dt z3 



0. 



longer than the interface amplitude Co- From Eq. ()2.6|) . we have V 
Upon neglecting an arbitrary function of time, which can absorbed into a redefinition 

of the velocity potential (j){x,z,t), we conclude that 



PPs d(j){x,z,t) , ^ ^ , „^ 

p(x,z,t) = ^ - Ps , do<z<d. (2.8) 

z-^ ot 



Note that the first term on the right side leads to an extremely large pressure for small 
z, which drives the formation of the solid "dead layer" of helium atoms close to the 
substrate. Incompressiblity of the superfluid gives us an additional equation for the 
velocity potential within the film, namely 

V20(x,z,t) =0. (2.9) 

Understanding third sound requires that we now specify the boundary conditions 
associated with Fig. 12.21 An undulating interface gives rise to a restoring force due to 
van der Waals interactions with the substrate in thin films, analogous to the earth's 
gravitational field for gravity waves in the ocean. Like oceanic gravity waves, it can 
be shown that the restoring force due to surface tension can be neglected at long 
wavelengths p4l [36]) . The velocity perpendicular to the wall of the inviscid fluid must 
vanish close to the dead layer, but some care is required with the upper boundary, 
which is itself a dynamical quantity, determined by the fluctuating interface position 
(i + C,{x,t). Upon assuming film undulations so slow that the pressure in the vapor 
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phase just above the film remains fixed at its equifibrium value po, we have from Eq. 
(EiHl) that 

Pp, d(t){x,d + C,t) 

TdTcr^'' — dt — =^°' ('-'^^ 

which upon linearization in C leads to ~^Jl^' — Ps ^'^^g^'^'*^ = Po — ■ The term on the 
right hand side can be eliminated by letting <f> (j) — {t/ps){po — /3ps/(P){this change 
of variables leaves Vg — unaffected), to obtain finally 

3/3 d(j){x,d,t) 

^C(a;, t) + g-^ = 0. (2.11) 



^ d(f>{x ,d,t) 

z=d+C 



A second equation associated with the free surface arises because the z-component of 
the velocity at the interfaces can be calculated as both v,- = ^'^^gl,^ 

and V2 = ^^g'^'*'' , which implies that 

d(p{x,d,t) d({x,t) 



dz dt 



(2.12) 



These two boundary conditions at the free surface can be combined into a single 
equation upon differentiating Eq. (|2.1ip with respect to time. We are left with a single 
set of equations for the velocity potential, namely 

V^4>{x,z,t) =0 (bulk) 

3p d4>{.,d,t) ^ d-^{x,d,t) (at the free surface) (2.13) 



(x,d 
"51 



^^^^^^^ - (near the substrate) 



This unusual wave equation with a parabolic time-dependent boundary condition 
has the solution 0(x, z, t) = (fiQ cosh(fcz)e'*''^~'"^'^^* in the limit kdo <C 1, with dispersion 
relation 

uj(k) = \^kta.nh(kd) ^ \^k = C3k, fcd<l. (2.14) 
V a \ d"^ 

The basic features of the "third sound" excitations with speed C3 described above have 
been well- verified in numerous experiments on helium films at low temperatures 
provided one replaces Eq. (|2.14p by 



w(fc)wy^^fc = C3fc, fcd<l, (2.15) 

where the factor Ps/p that enters the third sound veloicity C3 arises because only a 
fraction of the film is superfluid when the temperature T > (j27p . However, as the 
temperature gradually increases, a new source of dissipation damps out third sound 
waves, until they no longer propagate (|29l [30)) . 
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Superfluid vortices: analogy with 2d 
electrostatics 



The new source of dissipation for third sound/van der Waals waves in hchum films 
(which sets in well below the bulk superfluid transition temperature) is associated with 
a breakdown of the irrotational flow assumption embodied in Eq. (j2.5p . If we were 
analyzing gravity or capillary waves in a classical fluid, we could simply study the 
viscous dissipation associated with adding a continuous vorticity field to the potential 
flow discussed above. Vorticity can arise in superfluid helium films as well, but it is 
quantized. Indeed, the connection of the superfluid velocity with the phase of a single- 
valued quantum-mechanical wave function (see Eq. (j2.5p ) leads immediately to 

Vs{f) ■ dl^ — (h V6'(f) •dr= —27rs, s = 0,±1, ±2,... (3.1) 



If an array of vortices with topological "charges" {sj} occupy a set of points {rj}, this 
condition on the velocity field in an infinite plane leads to a nonzero 2d vorticity a;(r), 
namely 

uj{r) z- [V X Vs{t^)] ^ 27r—^SjS{f~fj). (3.2) 



For a single vortex with s = 1 located at tq = (xq, j/q) wc have d{x, y) = tan ^ ( ) , 
so that 

^.« = --|^,«."|^, (3.3) 

a velocity field that is curl- free except at r = rg . Insertion of this velocity field into 
the superfluid kinetic energy displayed in the second column of Table I, leads to a 
logarithmically diverging energy, E = 'Kps(h/m)'^\n-{L/a) , where L is the in-plane size 
of fllm, and a is a microscopic cutoff, of order the spacing between helium atoms. 
Henceforth, we restrict our attention to charge neutral vortex pairs, which have a 
flnite interaction energy. 

A detailed consideration of the damping of third sound by point-like vortex pairs in 
helium films (there is an interesting analogy with Maxwell's equations (HH)) would take 
us well beyond the scope of these lectures. To motivate our investigation of vortices 
on cylinders, we discuss here only the energetics of vortex pairs in the presence of a 
uniform supercurrent u^. Although a close approximation to this background superflow 
can be created via a torsional oscillator experiment (jSS]) . such a situation also arises 



Superfluid vortices: analogy with 2d electrostatics 11 



as the long-wavelength limit of the third sound waves discussed above: At wavelengths 
very much greater than the film thickness, the in-plane velocity field associated with 
the slowly-varying third sound wave looks approximately uniform with magnitude 
to a tightly-bound vortex pair, similar to, say, electron-hole pairs in a semiconductor 
exposed to electromagnetic radiation at optical wavelengths much larger than the pair 
separation ([1]). 

Energy calculations for helium films are facilitated by an analogy with 2d elec- 
trostatics, as first pointed out and exploited in the seminal work of Kosterlitz and 
Thouless p9)) . Kosterlitz and Thouless showed how defects such as dislocations and 
vortices could give rise to phase transitions in two-dimensional magnets, superfluids 
and crystals. Their work is directly applicable to the vortex-unbinding transition that 
causes the vanishing of the superfluidity density with a universal jump discontinu- 
ity in helium films (|40p at a two-dimensional critical temperature Tc, now called the 
Kosterlitz-Thouless transition. We will use similar ideas to determine how vortex pairs 
polarize in the presence of uniform supercurrents when < T ^ Tc- First define a 
fictitious electric field E{r) = zx Vs{x, y), which is simply a 90° counterclockwise rota- 
tion of the velocity field. Eq. ()3.2|) then becomes a two-dimensional version of Gauss' 
law, 

V ■ E{f) = 2T:—J2s]5{r-rj)- (3.4) 

j 

At long wavelengths, the coarse-grained 2d in-plane velocity field Vs{x,y), averaged 
over the film thickness, greatly exceeds the out-of-plane component, so that we now 
have a 2d incompressibility constraint, V • Ws = dxVg + dyV^ — 0. In terms of the 
fictitious electric field defined above, this constraint reads 

z-{VxE)=0, (3.5) 

which leads us to define an electrostatic potential $(a;, y) via E{x,y) = —W^{x,y) 
that satisfies 

V^^(x,y) = -2TT—y^ s,6(f-fj). (3.6) 
m ^-^ 

j 

In the continuum limit, ^(r) is related to the stream function of 2d fluid mechanics. 

As illustrated in Fig. 13. li this analogy with electrostatics maps computing the 
kinetic energy of the velocity configuration associated with a vortex pair immersed in 
a uniform supercurrent, which can be written as an integral over the intensity of the 
"electric field" 



1 /■ .o _ ,o 1 



2'^ / ' ' 2'^ m 



U^-Ps I d'r\v,\' = T^Ps— j d'r E{r) = ^Ps^ j d'r V<I>(f) , (3.7) 



2 



1 



2 



onto the 2d electrostatics problem of two charges in a uniform electric field proportional 
zxiP^. Let us define a "quantum of vorticity flux" = 2Trh/m. For two vortices located 
at fi and r2 with unit "charges" si = +1 and S2 = ^1, the electrostatic potential is 

$(f) = ^ ln(|f - fil/a) + ^ ln(|f - f^l/a) - {z x ««) • r. (3.8) 
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Fig. 3.1 (a) Vortex dipole with lines of constant pliase in a iielium film in the presence 
of a uniform vertical supercurrent (dashed lines), (b) The equivalent two-dimensional 
electrostatic problem consists of a pair of charges immersed a uniform horizontal electric field 
due to a charged capacitor, as the plate separation is taken to infinity. 

Upon inserting this potential function into Eq. ()3.7|) . we obtain, after integrations by 
parts and use of Eq. p.6p . the energy function for the vortex pair, 

^2 

U{ri - f2) = const. + p^-^ ln(|ri - r2\/a) - Ps^q{z x -D^) • (ri - r2), (3.9) 

where the constant depends on the detailed physics near the vortex cores. Although 
the logarithmic binding potential dominates when the vortex pair separation r = |fi — 
^2! ^ Tc = a;o/(27ru°). Fig. 13.21 shows that the linear contribution to the interaction 
potential J7(r) leads to a saddle point along the a;-axis at r = (re, 0). Beyond this saddle 
point, the pair can separate indefinitely. A detailed theory of this thermally activated 
"escape over a barrier" problem (pS)) (this theory now 1/3 of a century old!) reveals 
how long wavelength third sound waves are disrupted at low temperatures by this 
vortex unbinding process, with related results for the torsional oscillator experiments 
of Ref. ([38|). It is not our purpose here to recapitulate this ancient theory of vortex 
unbinding in planar helium films. Rather, we conclude these lectures by discussing 
how Eq. p.9p is altered, when we consider helium films on a cylinder. 
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Fig. 3.2 (Adapted from Ref. [26]) Contours of constant potential between a pair of oppositely 
charged vortices at separation r — {x,y). The potential decreases upon approaching origin, 
but has a short distance cutoff at radius a. A uniform superfluid velocity if^ has been imposed 
in the y direction. There is a saddle point at (rc,0), over which the pair separation can 
"escape" to large positive x. 
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Superfluid vortices on a cylinder 



Inspired in part by the interacting dislocation defects on cylinders discussed in the 
introduction, we now ask how quantized vortex pairs interact in the presence of a 
uniform supercurrent in a cylindrical geometry. What is the effect of the periodic 
boundary conditions, especially when the cylinder is long and narrow? A uniform 
supercurrent on a cylinder could be generated experimentally by first coating a solid 
spindle with a helium film, and then spinning up this cylinder at constant angular 
velocity at temperatures above superfluid phase transition, so that the fluid viscosity 
initially forces the film to follow the motion of the substrate. If the rotating spindle is 
subsequently cooled well below the transition, and the rotation stopped, the superfluid 
fraction of the film will continue to rotate indefinitely (|22p . The most important decay 
mechanism for this supercurrent will be nucleation and separation of vortex pairs, as 
illustrated in Fig. 14.11 Indeed, separation of a +/- vortex pair to infinity along the 
cylinder axis causes a 27r phase slip in the condensate order parameter one moves 
around the cylinder ((22l [26] |44)) , which reduces the background superfluid velocity, as 
discussed later in this section. 

We first treat the interaction of a vortex pair in the absence of a supercurrent. 
Suppose a positive vortex charge is glued to the origin, and its negative counterpart 
sits at position r= {x,y). We take the y-axis in the azimuthal direction, periodically 
repeated with period W = 27ri?o around the cylinder, and assume that the cylinder 
is infinite in the x-direction. We need to solve Eq. p.6|) . subject now to the periodic 
boundary conditions imposed by the cylinder. In the infinite plane, the energy would 
be C/(r) = Uo\n{r/a), where Uq = PsUJo/C^tt), with a corresponding attractive force 
exerted on the + vortex given by F = —Uof/r^. The generalization of this force to a 
cylinder requires summing an infinite number of negative image charges, 



(The forces due to the positive image charges cancel, as they must so that the positive 
charge at the origin does not exert a force on itself.) Similar sums over Matsubara 
frequencies are required to compute the Green's functions that arise quantum many 
body theory at finite temperatures (I4ip . and in many other fields. They can be carried 
out exactly using the Sommerfeld- Watson formula, which reads (jH]): 




(4.1) 




(4.2) 
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Fig. 4.1 Top: Vortex pair excitation in tlie presence of a constant background supercurrent 
ifg that wraps around tiie cylinder. Tliere is a force at riglit angles to this supercurrent (often 
called the Magnus force in fluid mechanics texts) that drives the vortices to separate towards 
the ends of the cylinder. The response to more general supercurrents that spiral around the 
cylinder is discussed in the text. Bottom right: periodic replicas of the positive and negative 
vortices required by the cylindrical boundary conditions. 

where the clockwise contour C for this complex contour integral encloses the entire 
real axis. The proof of this formula (valid for functions f{z) that fall off sufficiently 
rapidly for large z) proceeds by noting that the function cot(7rz) inside the integrand 
has only simple poles of residue unity which lie on the x axis at integer values. The 
sum is then evaluated by noting that since f(z) ~ l/|z| for large distances from the 
origin, we can deform the contour so that it captures only the poles of f{z) in the 
complex plane (note that the integral on the circle at infinity vanishes even though 
the decay is only ~ 1/|^|, due to the cot(7rz) term.) The results of this summation 
procedure for the forces are: 

F. = k^ {coth[7r(x - ty)/W] + coth[7r(a- + iy)/W]}, 

P _ ttUq sin(27ry/W) (4.dj 

~ W cos{2iry/W)-cosh{2Trx/W)- 

Note the periodicity in the y-coordinatc. Eq. (j4.3|) can be integrated to give the inter- 
action energy for two vortices on a cylinder, namely 

U{x, y) = const. + -Uq 

a result originally obtained by Machta and Guyer (j43|) . who modeled helium films 
in porous media as a network of interconnected cylinders. The contours of constant 
energy are displayed in Fig. 14.21 It is easily checked that Eq. (|4.4|) reduces to the 



W 

In i — sinh[7r(x + iy) /W] } + c.c 
na 



(4.4) 
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logarithmic interaction expected for the infinite plane when <^ W = 2t:Rq. 

More interesting, however, is the limit x ^ W, for which the forces in Eq. (|4.4p read: 

Fee ~ T^sgn(.T) + 2^e-2--/i^cos(2^2//M^), 



P ^ 27r;7( 



og-2WW'sin(27ry/PT/). ^ 



where sgn(x) is the sign of x. Note the similarity to the electric field between two 
parallel capacitor plates, with exponential corrections due to the discreteness of the 
charges. The cylindrical environment confines the electric field, so that it points pre- 
dominantly along the cylinder axis. Upon neglecting the exponential corrections and 
integrating, we see that the logarithmic binding potential for vortices is replaced by a 
linear binding for large vortex pair separations along the cylinder axis, 

U{x,y)^^\x\, \x\^W. (4.6) 

However, the potential function that tends to separate vortices in the presence of a 
circulating supercurrent like that in Fig. I4.1l is also linear in x\ (see Eq. p.9p ). Hence, 
it is far from clear at this point that the escape-over-a-barrier scenario summarized in 
Fig. 13.21 for an infinite planar film actually works for vortices on a cylinder, at least 
for small background supercurrents if^. Which linear potential wins at large x, the 
repulsive one due to the supercurrent, or the attractive one due to electrostatics of 
point charges on a cylinder? The remainder of this section provides a careful answer 
to this question, allowing for the more general case of a spiral supercurrent if^ inclined 
at a finite angle to the y-axis. 

The topological stability of supercurrents in the geometry of Fig. 14.11 can be un- 
derstood via the following argument: If i/jqIx,!/) is the single- valued condensate wave 
function describing to a superfiuid helium film at rest with respect to a cylinder, adi- 
abatically turning on a superfluid velocity Vs = {v^,vy) (which might in general be 
inclined at a finite angle to the y-axis) leads to a new wave function 

ip{x, y) = cxp[imu5 • r/h]^l:Q{x, y), (4.7) 

where f= (x, y)and m is the mass a helium atom. For notational simplicity, we now 
drop the superscript on the background superfluid velocity, and set if^ = Vg. The new 
wave function ip{x, y) must be single-valued as well; hence, we have ipi^-i V + ^) = 
'ip{x,y). The choice || x represents a supercurrent along the cylinder axis which 
is not quantized for an infinitely long cylinder. However, the choice \\ y describes 
an azimuthal supercurrent which must be quantized, according to u| = P — 

0,±1,±2, .... Hence, a general supercurrent must take the form 

vs=v:x+-'^py, p = 0,±l,±2,.... (4.8) 

Thus the allowed values in the (wf , wf )-plane are a set of horizontal lines with spacing 
^1^. As discussed in Refs. ((22|) and (|26p . vortex pair unbinding events mediate tran- 
sitions between various states of superfiow. A general supercurrent will spiral around 
the cylinder, leading to a set of complementary spiral trajectories for the optimal 
unbinding of vortex pairs. 
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Fig. 4.2 Equipotential contours for a positive vortex witii unit cliarge on a cylinder, inter- 
acting with a negative vortex at the origin. The interaction is logarithmic for x <^ 2hRq. 
However, for x ^ 2nRo, the energy contours are parallel to the y axis and reflect a linear, as 
opposed to a logarithmic, binding potential. 

Is there, similar to Fig. l3.2l in the infinite plane, always a saddle point for the escape- 
over-a-barrier energy landscape representing a vortex pair superimposed on a nonzero, 
uniform spiral background supcrcurrent? Or, as one might naively guess, does a set of 
very small but nonzero supercurrents exist such that that the repulsive potential due 
to the background "electric field" cannot overcome the electrostatic linear potential 
for a vortex pair that tries to separate along the cylinder and the vortex pair remains 
forever bound? To answer this question, we need to understand the cylindrical energy 
landscape that generalizes Eq. ()3.9p . namely 

C/(x,y) = ic/o 

subject to the quantization condition (|4.8p above. We have suppressed a short distance 
constant contribution that can be incorporated into the vortex core energies. It is 
straightforward to analyze two limiting cases: 

Case 1 : t/g |j i (no quantization of on the infinite cylinder). We have z • (wg x r) = 
VsV and assume a saddle point of the form (0,j/c), found by minimizing U{0,y) = 
C/q ln[^ sin(7r?//iy)] — psio^VsV- It is easy to see that yc is the solution of cot(-^) = 
^VsW, which exists for any value of Vs- The solution ydvs) has the following prop- 
erties: (a) as Vs — >■ c», j/c 0, so that vortex unbinding is very easy; (b) as W ^ oo, 
we recover the infinite plane result, yc — rc ~ h/mvs — uiq/2t:Vs] and finally (c) 




W 

In { — sinh[7r(.T + iy)/W] 
Tia 



+ c.c 



(4.9) 
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lim rjc = W/2. This last result is to be expected, since vortex pairs can always unbind 

-Us— 

by going azimuthally around the cylinder, due to the periodic boundary conditions. 
Thus, as we let f ^ — !• along the x-axis (continuously, because there is no quantiza- 
tion), we always have a saddle point, and the conventional picture of Refs. ([221) and 
(j26p for the gradually destruction of superfluidity by successive unbindings of vortex 
pairs will hold. 

Case 2 : Vs || y. We now have a quantized superfluid velocity: = -^^P, P = 
0, ±1, ±2, .... In this case z ■ (w, x f) = —VsX and we assume a saddle point of the form 
(xc,0). Consider a directed along the negative y-axis, which requires that wc now 
minimize U (0, y) — Uq ln[-^ sinh(7ra;/W^)] — pstooVsX. Although there is no saddle point 
for p = 0, there is in fact a saddle point for every nonzero value of p, given by the 

solution of coth(^^) = 2p, or explicitly by Xc{p) = ^ In (^ 2p-i ) • -Note in particular 

that x™'" = Xc{p)\p=i = l^lnS. For large W we have, lim Xc{p) = = h/mvs ~ 

cjo/SttWs, as in flat space . 

The "nai've" conclusion that no saddle point exists for small but nonzero supercur- 
rents in the y-direction on the cylinder arises from comparing two competing linear 
terms in the vortex pair potential. 



lim ;7(x,0) 



— -ps^OVs 



(4.10) 



which leads to a "nai've" critical current v': = = i — Itt. However, this is half 

the minimum allowed nonzero azimuthal supercurrent on a cylinder, and hence cannot 
not be realized experimentally, due to the constraint of a single- valued wave function. 

The existence of a saddle point for cases 1 and 2 above strongly suggests that a 
saddle point exists for any choice of Vs- This is indeed the case, as we shall now prove: 
Since the potential energy of the vortices has an exact electrostatic analogy, it is clear 
that the potential energy is a harmonic function, and thus any extremum has to be 
a saddle point. Let us assume without loss of generality thatu^ > 0,v^ < 0. Using 
the arguments outlined above, we know that for x ^ W and fixed y the potential is 
approximately linear in x, and the slope must be negative since the argument relying 
on the quantization of still holds. Consider now a set of arbitrarily chosen paths 
(not to be confused with the energy contours or the unbinding trajectories!) each of 
which start at the origin and ends at a; — >■ oo, as shown in Fig. 14.31 The potential 
energy near the origin is singular and attractive, so the slope as one moves from the 
origin towards positive x is positive. Thus, each of the paths has to have a maximum 
somewhere along it. Moreover, from the properties of the potential we shall now show 
that each path must have a single maximum along it: consider the x component of the 
force that a vortex would encounter were it to be dragged along one of these paths. 
For any extremum this component must vanish. Along the path, the force component 
due to the supercurrent is constant, while the force due to the interaction with the 
vortex at the origin is monotonically decreasing (since it can be thought of as the 
superposition of the vortex and its image charges, all of which are monotonically 
decreasing). Therefore there cannot be more than a single maximum. 
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Fig. 4.3 A dense set of paths on the cyhnder. The figure is an illustration used in the proof 
of the existence of a saddle point for any direction of the supercurrent. For each path the 
maximum of the potential energy is chosen, resulting in the line joining points A and B. It 
can now be proved that there is a minimum on this curve, which is a saddle-point of the 
potential energy. 

For the lines x=0 and y=0, these maxima (points A and B in the figure) were 
previously discussed. Furthermore, wc know that at point A the force due to the 
vortex at the origin has no component in the y direction (from symmetry), such that 
the gradient of the total potential energy is only contributed by the supercurrent and 
therefore points in the positive y direction. Similarly, the gradient at point B is in the 
positive X direction. These are denoted by dashed arrows in Fig. 14.31 From continuity 
of the potential, and using the fact that there is a single maximum for each path, the 
set of maxima for all of the paths must lie on a single continuous path, denoted in the 
figure by the curve joining points A and B. Now let us consider the minimum of the 
potential along this curve. It cannot lie on the boundary, since the gradients at the 
edges point inwards. Therefore it must lie at some interior point. This point is thus 
a minimum in one direction, and a maximum in another - and is therefore a saddle 
point. 



5 

Dislocations on cylinders are 
different 



As discussed in Sec. I, dislocation defects in two dimensional elastic solids, whose close 
analogy with vortices in helium films is highlighted in table I, play an important role 
in in the elongation of the cell walls of cylindrical bacteria. Fig. 11.21 summarizes the 
model of bacterial elongation explored in Refs. ([T5)) and Another example of 2d 
crystals with defects on a cylinder arises for interacting colloids on the surface of a 
liquid film coating a solid cylinder, where repulsive short-range forces give rise to the 
self-organized emergence of a two dimensional crystalline solid. Defects in colloidal as- 
semblies on the related curved surfaces of capillary bridges have been recently studied 
experimentally and theoretically (j45p . Here, the Gaussian curvature can be positive or 
negative; the zero Gaussian curvature of a cylinder is a special case Fig. lS.ll shows 
an alternative dislocation pair, with equal and opposite Burgers vectors ±b now par- 
allel to the azimuthal direction. This configuration might allow a bacterium to thread 
through a narrow circular pore smaller than its natural ^1/xm diameter, to which it 
relaxes on the left and right sides of the pore in the figure. The circular constriction 
produces a two-dimensional stress (force per unit length) ayy that tends to drive the 
dislocations apart. As discussed below, the dislocations are logarithmically bound for 
distances short compared to the circumference of the cylinder, similar to a ± vortex 
pair. The indicated component ayy of the stress tensor tends to squeeze the horizontal 
lattice rows together. Like the vortex pair in the presence of a background supercur- 
rent shown in Fig. 14.11 this force leads to an additional linear potential forcing the 
two dislocations apart. Once these defects separate completely (for this geometry, an 
entire row of unit cells would have to be removed along the cylinder) , the crystal would 
have reduced its circumference by one lattice constant. Superficially, the dislocation 
configurations such as those shown in Figs. 11.21 and 15.11 seem quite similar to super- 
fluid vortices on a cylinder. Nevertheless, as highlighted below, there are important 
and surprising differences, due to the vector nature of the Burgers vector topological 
charges of the dislocations. See Ref. ([20| for a more thorough and extensive discussion 
of the physics of dislocations on cylinders summarized briefly below. 

In the infinite plane, the interaction energy of two edge dislocations with Burgers 
vectors bi and 62, and with a relative separation r = fi — r2, in the presence of a 
uniform background stress tensor aij , is given by 

- bhCTjkSijzn, (5.1) 



U{r) = const. + 



Y 
An 



bi-b2H\r\/a) 



(6i-r)(62-r) 
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Fig. 5.1 (adapted from Ref. [18]) Illustration of how dislocations could facilitate remodeling 
of a cylindrical bacterium growing through a narrow constriction bounded by the vertical 
dashed line segments. The cell wall assumes its natural radius i?o of about 1/xm outside the 
central region, which exerts a stress Oyy on the bacterium in the azimuthal direction. If the 
dislocations are able to separate in the horizontal direction, the radius of the cell wall will 
shrink to accommodate the pore. 

where a is a microscopic cutoff of order the lattice spacing and Y is an elastic constant, 
the 2d Young's modulus for the case of an isotropic crystal ([7]). The last term is the 
energy corresponding to the Peach-Koehler force ((SJ HTj) due to the stress field ajk , 
where Sijz is the Levi-Civita tensor. Comparison with Eq. p.9|) for the energy of 
a vortex pair with a background supercurrent shows that for crystalline solids, the 
Peach-Kohler force is the analog of the Magnus force on a vortex. In a coordinate 
system where aij is diagonal {axy = 0), this force is -F = —haxxV when b ~ bx (as in 
Fig. II. 2p and F = bcjyyX when b = by (as in Fig. 15. 

On a cylinder, the potential energy of a dislocation pair takes the form 

U{r) = const. + E{bi,b2; x,y) - bkcrjk£i]zri, (5.2) 

where the linear potential due to a background stress is unchanged, and the inter- 
action energy E{bi,b2; x,y) is periodic in y with period W = 2'kRq and can be cal- 
culated by summing over periodic image dislocations using the Sommerfeld- Watson 
method as sketched above for vortices ((20|) . Fig. 15.2b , shows the energy contours for 
E{by,—by;x,y), i.e., for dislocations oriented as in Fig. 15.11 The potential is isotropic 
and logarithmic for short distances, and the dislocation pair experiences a linear po- 
tential at large separations along the x-axis, like vortices on a cylinder. However, the 
energy contours for E(bx,—bx;x,y), i.e., for dislocations oriented as in Fig. 11.21 dif- 
fer dramatically, as shown in Fig. 15.2b . The pair interaction, although isotropic and 
logarithmic for small separations, is now exponential screened when x ^ W, i.e. for 
separations much larger than the circumference lf20| . The source of this remarkable 
behavior is the dipole-like angular interaction term in Eq. (|5.1[) on the cylinder. When 
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Fig. 5.2 (a) Equipotential contours for an interacting dislocation pair, one at the origin 
and another at {x,y), on a cylinder with azimuthal Burgers vectors similar to Fig. 15.11 (b) 
Similar plot for a dislocation pair with Burgers vectors parallel to the cylinder axis, as in 
Fig. 11.21 Unlike the linear potential at large horizontal separations in (a), the interactions are 
exponentially screened in this case. 

b is parallel to x, as is the case for the dislocations that mediate bacterial elongation in 
Fig. ll.2[ the angular term causes the long range strain field from an isolated dislocation 
to fall off exponentially fast. The reason for this strange behavior is that isolated dislo- 
cations on a cylinder with b\\x act like grain boundaries in an infinite crystalline solid, 
due to the periodic boundary conditions. Grain boundaries connect two orientationally 
mismatched crystals, and can be regarded as an infinite row of aligned dislocations 
with the same topological charge ([5l [71 148)) . The innocuous- looking angular term in Eq. 
(jS.ip leads to a low energy cost for grain boundaries, contrary to what one might have 
guessed by a nai've application of the electrostatic analogy to a array of dislocations 
with identical charges. The vector nature of these charges is crucial. 

Thus, despite the many deep and useful analogies between superfluid helium films 
and 2d crystalline solids, important differences emerge, due to the vector nature of the 
topological charges characterizing the dislocations. See Ref. pO)) for more details. 
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